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Abstract. In this paper, we investigate frames for L^— 7T, ir] d consisting of exponential 
functions in connection to oversampling and nonuniform sampling of bandlimited func- 
tions. We derive a multidimensional nonuniform oversampling formula for bandlimited 
functions with a fairly general frequency domain. The stability of said formula under 
various perturbations in the sampled data is investigated, and a computationally man- 
agable simplification of the main oversampling theorem is given. Also, a generalization of 
Kadec's 1/4 Theorem to higher dimensions is considered. Finally, the developed techniques 
are used to approximate biorthogonal functions of particular exponential Riesz bases for 
L>2[— 7r, 7r], and a well known theorem of Levinson is recovered as a corollary. 

1. Introduction 

The subject of recovery of bandlimited signals from discrete data has its origins in the 
Whittaker-KotePnikov-Shannon (WKS) sampling theorem (stated below), historically the 
first and simplest such recovery formula. Without loss of generality the formula recovers a 
function with a frequency band of [—it, tt] given the function's values at the integers. The 
WKS theorem has drawbacks. Foremost, the recovery formula does not converge given 
certain types of error in the sampled data, as Daubechies and DeVore mention in [7J. They 
use oversampling to derive an alternative recovery formula which does not have this de- 
fect. Additionally for the WKS theorem, the data nodes have to be equally spaced, and 
nonuniform sampling nodes are not allowed. As discussed in [151 pages 41-42], nonuniform 
sampling of bandlimited functions has its roots in the work of Paley, Wiener, and Levin- 
son. Their sampling formulae recover a function from nodes (t n ) n , where (e ltnX ) n forms a 
Riesz basis for Z^f— 7r,7r]. More generally, frames have been applied to nonuniform sam- 
pling, particularly in the work of Benedetto and Heller in [2] and [3]; see also |15|. chapter 10]. 

In Section 3, we derive a multidimensional oversampling formula, (see equation (4)), for 
nonuniform nodes and bandlimited functions with a fairly general frequency domain; Sec- 
tion 4 investigates the stability of equation (4) under perturbation of the sampled data. 
Section 5 presents a computationally feasible version of equation (4) in the case where the 
nodes are asymptotically uniformly distributed. Kadec's theorem gives a criterion for the 
nodes (t n ) n so that (e LtnX ) n forms a Riesz basis for L,2[— vr, tt]. Generalizations of Kadec's 1/4 
theorem to higher dimensions are considered in Section 6, and an asymptotic equivalence 
of two generalizations is given. Section 7 investigates approximation of the biorthogonal 
functionals of Riesz bases. Additionally, we give a simple proof of a theorem of Levinson. 



i 



2 BENJAMIN BAILEY 

This paper forms a portion of the author's doctoral thesis, which is being prepared at Texas 
A & M University under the direction of Thomas Schlumprecht and N. Sivakumar. 

2. Preliminaries 
We use the d-dimensional L2 Fourier transform 

Hf)(-)= [ /(£)< f£L 2 (R) d , 
where the inverse transform is given by 

(2vr) a! J Rd 

This is an abuse of notation. The integral is actually a principal value where the limit is in 
the L2 sense. This map is an onto isomorphism from L2(M. d ) to itself. 

Definition 2.1. Given a bounded measurable set E with positive measure, we define 
PW E ■= {/ G L 2 (M d )|supp(7'~ 1 (/)) C E 1 }. Functions in PW E are said to be bandlim- 
ited. 

Definition 2.2. The function sine : R — > R is defined by sinc(x) = sin ^ . We also define 
the multidimensional sine function SINC : R d -»• M d by SINC(s) = sinc(xi) • . . . • sinc(cc d ), 
ac = (aci, . . . ,x d ). 

We recall some basic facts about PWe- 

1) PWe is a Hilbert space consisting of entire functions, though in this paper we only regard 
the functions as having real arguments. 

2) In PWe, L2 convergence implies uniform convergence. This is an easy consequence of 
the Cauchy-Schwarz inequality. 

3) The function sinc(7r(x — y))) is a reproducing kernel for PWr_ f j, that is, if / G PW^_ n ^, 
then we have 

/oo 
/(r)sinc(vr(t - r))dr, t G R. 
-00 

4) The WKS sampling theorem (see for example [14., page 91]): If / G PWr^ j, then 

/(*) = ^/(n)sinc(vr(t-n)), tei, 

where the sum converges in PW^_ W ^, and hence uniformly. 

If (/n)neN is a Schauder basis for a Hilbert space ff, then there exists a unique set of func- 
tions (fn)n<=N, (the biorthogonals of (/ n ) n eN) such that (f n , fm) = ^«m- The biorthogonals 
also form a Schauder basis for H. Note that biorthogonality is preserved under a unitary 
transformation. 



SAMPLING AND RECOVERY OF MULTIDIMENSIONAL BANDLIMITED FUNCTIONS VIA FRAMES 3 

Definition 2.3. A sequence (f n ) n C H such that the map Le n = f n is an onto isomorphism 
is called a Riesz basis for H. 

The following definitions and facts concerning frames are found in [BJ section 4]. 

Definition 2.4. A frame for a separable Hilbert space H is a sequence (/ n )n C H such 
that for some < A < B, 

(2) A||/|| 2 <^|(/,/„)| 2 <i?||/|| 2 , V/etf. 

n 

The numbers A and B in the equation ([2]) are called the lower and upper frame bounds. 

Let H be a Hilbert space with orthonormal basis (e n ) n . The following conditions are equiv- 
alent to (/n)n C iT being a frame for f/\ 

1) The map L : H —■ H defined by Le n = f n is bounded linear and onto. This map is called 
the preframe operator. 

2) The map L* : H — > H (the adjoint of the preframe operator) given by / h- > ^2 n (f, fn)en 
is an isomorphic embedding. 

Given a frame (f n ) n with preframe operator L, the map S = LL* given by Sf = ^2 n (f, f n )fn 
is an onto isomorphism. S is called the frame operator associated to the frame. It follows 
that S is positive and self-adjoint. 

The basic connection between frames and sampling theory of bandlimited functions (more 
generally in a reproducing kernel Hilbert space) is straightforward. If (e** n ^) n is a frame 
for / 6 -FW[_,r )7r ] with frame operator 5, and / € FWI-^j, then 

S(J r ~ 1 (f)) = '^2{J r ~ 1 (f), fn)fn = Y,H^\f)){tn)fn 

n n n 

implying that J r_1 (/) = En /(*n)S _1 /n, so that / = £ n f{t n )HS~ X fn)- Note that in the 
case when i n = n, we recover the WKS theorem. 

Definition 2.5. A sequence (/ n )n satisfying the second inequality in equation ([2]) is called 
a Bessel sequence. 

Definition 2.6. An exact frame is a frame which ceases to be one if any of its elements is 
removed. 

It can be shown that the notions of Riesz bases, exact frames, and unconditional Schauder 
bases coincide. 

Definition 2.7. A subset S of W 1 is said to be uniformly separated if 

inf \\x — t/ | j 2 > 0. 

Definition 2.8. If S = (xk)k is a sequence of real numbers and / is a function with S in 
its domain, then fs denotes the sequence (/(x/c))fc- 
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3. The multidimensional oversampling theorem 
In [7j, Daubechies and DeVore derive the following formula: 

(3) /(^^^-l)' tGM ' 

where g is infinitely smooth and decays rapidly. Thus oversampling allows the representa- 
tion of bandlimited functions as combinations of integer translates of g rather than the sine 
function. In this sense equation (|3j) is a generalization of the WKS theorem. The rapid 
decay of g yields a certain stability in the recovery formula, given bounded perturbations 
in the sampled data [7]. 

In this section we derive a multidimensional version of equation ([3]), (Theorem 3.1) for un- 
equally spaced sample points, and the corresponding non-oversampling version of the WKS 
theorem is given in Theorem 3.2. 

Daubechies and DeVore regard J-~ l (f) as an element of L2\ — Xtt, Xtt] for some A > 1. In 
their proof the obvious fact that [— it, it] C [—Xtt, Xtt] allows for the construction of the bump 
function J 7 " 1 (3) G C°°(R) which is 1 on [— vr, 7r] and off [— Xtt, Xtt]. If their result is to be 
generalized to a sampling theorem for PWe in higher dimensions, a suitable condition for 
E allowing the existence of a bump function is necessary. If E C R rf is chosen to be compact 
such that for all A > 1, E C int(A-E'), then Lemma 8.18 in page 245], a C°°-version of 
the Urysohn lemma, implies the existence of a smooth bump function which is 1 on E and 
off XE. It is to such regions that we generalize equation ([3]): 

Theorem 3.1. Let G E C R d be compact such that for all X > 1, E C int(AE'). Choose 
S = (^n)neN C R d such that (/ n )neN; defined by f n {- ) = e l (' ,tn ' , is a frame for L^iE) with 
frame operator S. Let A > 1 with J"" 1 ^) : R d -> R, T^ig) G C°° where T~ l (g)\ E = 1 
and J r ~ 1 (g)\(\ o E) c= 0- Lf X > Xo and f G PW E , then 

(4) /w = pE(E^/(t))^-|)> t£< 

fceN neN 

where B^n = (S~ l f n ,S~ 1 fk} E - Convergence of the sum is in L2(M. d ), hence also uniform. 
Further, the map B : £ 2 (N) -4 £ 2 (N) <ie/ined by (y k )keN ^ ( E«eN B k n yn) km is bounded 
linear, and is an onto isomorphism iff (f n )neN * s a Riesz basis for L2{E). 

Proof. Define f\ in (- ) = fn(j)- Note that (fx,n)n is a frame for L^XE) with frame operator 
Step 1: We show that 

(5) f = Y,f(^)^ s x l h,n)F-\a% f^PW E . 

n 
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We know supp(J" 1 (/)) C E C XE, so we may work with T 1 (/) via its frame decomposi- 
tion. We have 

T- 1 (f) = S^ 1 S x (T- 1 (f)) = Y,(F- 1 (f),fx,n)xES^ 1 fx, n , on XE. 

n 

This yields 

^" 1 (/)=E(- 7r_1 (/),/A 1 n)A E (^ 1 /A,n)^- 1 (ff), on R d , 
n 

since suppj r (g) C XE. Taking Fourier transforms we obtain 

(6) / = ^(^ 1 (/),/A,n)A^[(^ 1 /A,n)^- 1 (5)], On R d . 

n 

Now 

(F- 1 (f),hn)xE= [ ^- 1 (/)(e)e-^^^ = /(^ 
JXE v A 

which, when substituted into equation ([6]), yields ©. 
Step 2: We show that 

(7) f(-) = Y,f{j)[E^ 1 hn,S- x 1 h, k ) XE g(-- t -± 

n k 

where convergence is in L<i- 

We compute TKS^ 1 fx tn )J r ~ 1 (g)]. For h € L2(XE) we have 

ft, = S\(S X 1 h) = ^(<S* A 1 h, fx,k)\Efx,k = y^(fo; ^ 1 fx,k)xEfx,k- 
k k 

Letting /i = S^Va,™ 

S^fx^n = /^jSx 1 fx,n, 1 fx,k)XEfx,k- 
k 

This gives 

FKS^fx^F-Hg)^) = Y,( S x 1 f^ S x 1 f^)xEHfx,kJ'- 1 (g)](-) 



Y,(S^h,n,Sx- l fx,k)xE [ e^T-'igme-^dt 
^2(S^fx, n ,S^fx, k )xE [ J-^XOe^-x^ 

u JXE 



^A^X 1 /A,n ) S\ 1 fx,k)xEg{ r-), 



so ([7]) follows from ([5]). 
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Step 3: We show that 

(8) (S^fx,n,S^ 1 fx,k)xE = ^(S- 1 fn,S- 1 f k )E, for n, k G N. 

First we show (S^Va^X - ) = jd (S^VnX x)> or equivalently that / A ,„ = ^^((-S^/nXi)) ■ 

We have for any g € L%(XE), 

(9,fx,k)xE = [ g(Oe- i ^dti = \ d [ g(Xx)e-^dx = X d (g(X(-)),f k ) E . 

JXE JE 

By definition of the frame operator S\, 

S\9 = E^' fx,k)xEfx,ki 
fceN 

which then becomes 

Sxg = x d Y,(g(H-))Jk) E fx,k- 

k 

Substituting g = ^-(S 1-1 f n )(j) into the equation above we obtain 

^ A ((5" 1 /n)(^)) = Y,(S- 1 f n ,f k ) E fx,k= (^(S-Vn))^) =fx,n- 

k 

We now compute the desired inner product: 

(S^fx,n,S^f x ,k)xE = ^3 JjS^f^iS-ifk^dx 

= ^j E (S- 1 f n )(x)(S^f k )(x)dx = ^(S- 1 f n ,S- 1 f k ) E . 
Note that equation ([7]) becomes 

0) /(• ) = y d E f(j) [ E^A, -|)] . 

n k 

Step 4: The map V : £ 2 (N) 4(N) given by x = (x k ) keN i-> ( £) n B kn x n ) km = Bx is 
bounded linear and self- adjoint. 

Let (d k ) ke fq be the standard basis for ^(N), and let (e k ) ke fq be an orthonormal basis for 
L 2 (E). Then 

Vdj = (B kj ) keN = ^Bkjdk = '^2(S~ 1 f j ,S~ 1 f k )d k = ^(L* (S^ 1 ) 2 Lej , e k )d k , 

k k k 

where L is the preframe operator, i.e., S = LL*. Define <f> '■ ^ 2 (N) — > L 2 {E) by (p(d k ) = e k , 
k G N. Clearly is unitary. It follows that V = L* (S~ 1 ) 2 L(f>, which concludes Step 4. 
From here on we identify V with B. Clearly B is an onto isomorphism iff L and L* are 
both onto, i.e., iff the map Lc n — f n is cin onto isomorphism. 
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Step 5: Verification of equation (pp. Recalling Definition I2.8| f s / x = (/(x)) n eN' ^ or eacn 
t € W 1 , let g x (t) = (ff(*-^)) neN - Noting that G L 2^E), and recalling that 

{fxn)n is a frame for L 2 (A£'), we have 



(10) 
and 



£ M* " j) I' = E K^M* - j)),fx,n)x E ? < A x \\T-\g{t - -))| 



Note that equation ([9]) becomes 

/(*) = 4 £/(t) 



A rf ^ ' v A 

n fc 



■^E/(j)[E«-»(«-j) 

= ^(fs/x,BgMt)) = ^(Bf s/X ,g^{t)) 

n 



!) = ^E(Ew(^)>(«-|), 

fc6N neN 



which proves Q. 

Step 6: We verify that convergence in equation @ is in L20R) (hence uniform). Define 

= 4 E (J3/5/A)^(t-y) 



Kk<n 



and 



Then 



fm,n(t) = ^ E ( B /SAM*-T)- 



m<k<n 



[r-Hfm,n)](0 = yjE (S/fi/Ajfc^bC-y)] 



m<k<n 



m<k<n 



SO 



[^" 1 (/m J n)]||l = 4 / ^-"(^(Ol 2 ! E (^/^/a)^' 
A •/ A.E I rr ^>„ 

- Xd| E ( B fs/x)kf. 



m<k<n 
2 
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If (h n ) n is a orthonormal basis for L,2{\E), then the map Th k = f\k (the preframe operator) 
is bounded linear, so 

Wl?' 1 (f m ,n)]§ < ^\\t( £ (Zfs/xhhk) '<^I|T|| 2 Y, \( B fs/x)k\ 2 - 



m<k<n m<k<n 



But Bfs/x G ^ 2 (N), so ([[J 7 1 (/ m ,n)]||2 —> as m, n — > oo. As J 7 1 is an onto isomorphism, 
we have ||/ m ,n|! — >• 0, implying that ||/ — / n || — > as n — > oo. □ 



Note that equation (|3.ip is conveniently written as 

(11) f(t) = y d Y(Bf s/x ) k g(t- t f), teR d . 

k 

Remark: There is a geometric characterization of sets Bcl (i such that E C vak{XE) for 
all A > 0. Intuitively, E must be a "continuous radial stretching of the closed unit ball". 
This is precisely formulated in the following proposition (whose proof is omitted). 

Proposition 3.2. If G E cM. d is compact, then the following are equivalent: 
1) EC int(AE) for all A > 1. 

i2j There exists a continuous map 4> : S^ 1 -> (0,oo) such that E = {ty<p(y)\y G S d -\t G 
[0,1]}- 

The following is a simplified version of Theorem 13.11 which is proven in a similiar fashion: 

Theorem 3.3. Choose (t n )neN C M d such that (/ n )neN> defined by /„(■ ) = , w 2 e ^'' tn ^ * s 
a frame for Z^d - T) 7r] d ). If f £ PWe, then 

(12) /(t) = ^(^)S fcn /(t n ))siNC(7r(t-t fc )), 

fceN neN 

XTze matrix B and the convergence of the sum are as in Theorem \3.1\ 



Equation Q generalizes equation (TT2]) in the same way that equation © generalizes the 
WKS equation. 

We can write equation (fT2|) as 

(13) f(t) = ^2(Bf s ) k SWC(ir(t - t k )). 

keN 

The preceding result is similar in spirit to Theorem 1.9 in [H page 19]. 

Frames for L2(E) satisfying the conditions in Theorems 13.11 and 13.31 occur in abundance. 
The following result is due to Beurling in [5j see Theorem 1, Theorem 2, and (38)]. 

Theorem 3.4. Let A C R. d be countable such that 

r(A) := - inf ||A - /u|| 2 > 
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and -R(A) := sup inf ||A — /x[ [2 < tt- 

If E is a subset of the closed unit ball in M d and E has positive measure, then {e^'' A ^ |A £ A} 
is a frame for L2 (E) . 



4. Remarks regarding the stability of Theorem 13.11 

A desirable trait in a recovery formula is stability given error in the sampled data. Suppose 
we have sample values f n = f(j) + e n where sup n |e n | = e. If in equation ([3]) we replace 
/(t) by f n , and call the resulting expression /, then we have 

\f(t) - /(f) | < e i ^ \g(t - ~)| < eiX-'Wg'h, + \\g\W). 

It follows that equation ([3]) is certainly stable under perturbations in the data, while 
the WKS sampling Theorem is not. For a more detailed discussion see |7J. 

Such a stability result is not immediately forthcoming for equation as the following 
example illustrates. 

Restricting to d = 1, let (tn)nez satisfy to = D ^ Z, and t n = n for n ^ 0. The forthcoming 
discussion in Section 5 shows that (f n )nez is a Riesz basis for L2[— 7r,7r]. 

Note that when (f n ) n is a Riesz basis, the sequence f n ) n is its biorthogonal sequence. 
We matrix B associated to this basis is computed as follows. 

The biorthogonal functions (G n ) n6 z for (sinc(7r(- — n))) ne z are 

(n — D)[t — n) 
G (t) - S[Uc{7Tt) 



sinc(7rD) 

That these functions are in PW^j is verified by applying the Paley- Wiener Theorem 
|144 page 85], and the biorthogonality condition is verified by applying equation ([1]). Again 
using equation (P), we obtain 

D(-l) m 

i) B m0 = (G 0l G m ) = . - K ' — , m + 0, 

smc(7riJj(m — u) 

ii) B 00 = (Go, Go) ' 



sinc 2 (-7r.D) ' 



Hi) B mn — (G n , G m ) — Onm ~i~ 7 rvTT F>T> else. 

(n — D)(m — D) 
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Note that the rows of B are not in t\ , so that as an operator acting on , B does not act 
boundedly. Consequently, the equation 

(14) /(*) = j£(*/s/aM*-t) 

k 

is not defined for all perturbed sequences f$/\ where (fs/x)n = (fs/\)n+£n where sup n \e n \ = 
e. 

Despite the above failure, the following shows that there is some advantage of equation @ 
over equation (fl2jh 

If fs/\ is some perturbation of fg/\ such that \\Bf s j\ — Bf s /\\\oo < e, then 




5. Restriction of the sampling Theorem to the case where the exponential 

frame is a rlesz basis 

From here on, we focus on the case where (i n ) n eN is an perturbation of the lattice Z d , 
and (/ n )neN is a Riesz basis for L2[— it, 7r] d . In this case, under the additional constraint 
that the sample nodes are asymptotically the integer lattice, the following theorem gives a 
computationally feasible version of equation @ • The summands in equation Q involves 
an infinite invertible matrix B, though under the constraints mentioned above, we show 
that B can be replaced by a related finite-rank operator which can be computed concretely. 
Precisely, one has the following. 



Theorem 5.1. Let (n k ) k ^ be an enumeration ofL d , and S = (ifc)fceN C such that 



lim \\n k - ifclloo = 0. 

k— voo 



Define e k , f k : R d -> C by e k (x) = j^y^e^^ and ^372 e^*'^, and let (h k ) k be 
the standard basis for ^(N). Let Pi : ^(N) — > £2^) be the orthogonal projection onto 
span{/ii,-- - ,hi}. //(/fc)fceN is a Riesz basis for Z^f— 7r,7r] d , then for all f £ FlVr^jj, we 
have 



(16) f(t)= lim ^Y.^^'^fs/xhgit-^), t€R d , 

k=l 

where convergence is in L2 and uniform. Furthermore, 

(pu-ip\ _ / sinc7r(t ni i -t m>1 ) ■ ... ■smcir(t nid -t m>d ), l<n,m<l 
W l)nm ~ \ 0, otherwise. 

Convergence of the sum is in L2 and also uniform. 
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There is a slight abuse of notation in the formula above. The matrix PiB~ 1 Pi is clearly not 
invertible as an operator on I2, and it should be interpreted as the inverse of an I x I matrix 
acting on the first I coordinates of fs/x- 

The following version of Theorem 15.11 avoids oversampling. Its proof is similar to that of 
Theorem 15.11 



Theorem 5.2. Under the hypotheses of Theorem \5.1l 

1 

(17) f{t) = lim ytafl-^rVslfcSENTCCt _ tk ), t G R d , 

k=l 

where convergence of the sum is both L2 and uniform. 

The following lemma forms the basis of the proof of the preceding theorems, as well as the 
other results in the paper. 

Lemma 5.3. Let (n k ) k& ^ be an enumeration of T, d , and let (ijt)fcgN C M. d . Define e k , fk '■ 
R d -> C by e k (x) = i^e*^*'^ and f k (x) = jA^e^^K Then for any r,s > 1, and 

any finite sequence (a k ) s k=r , we have 
(18) 



ak j((-),n k ) _ _2h_ P m,t k ) 

(2vr)<V 2 (27r) rf / 2 



< f e ^r<k<J n k-tk\\oo) -{\(^T 



k=r 

Proof. Let 5 k = t k - n k where 5 k = (Ski, ■■ , o~kd)- Then 



a k 

k=r 



fc=r ^ ' k=r ' 

Now for any 5 k , 

31=0 JV j d =0 3d ' 

ji! ■ ■■■ ■ Jd! 

Ji>0 



.■.7i+-+.7„ (toi)-* 1 • ■ ■ ■ • (£fc^rf) Jd 

/ 

where J = {(ji, • • • ,jd) £ > 0, (ji, • • • ,jd) 7^ 0}. Then equation (fT9j) becomes 



M*> = -Et^^'I E ' 



h+-+3d ( 6 klXl) n ' ■ ■ ■ • (5 k dXd) 3d 



h } - ■ ■ ■ ■ -jd) 
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iA ■ 

(h,-,jd)£J 



Ji'- ■ ■■■ -Jd ] - 



/ j /o \rllo k\ ' ' ' kd 



k=r 



(2vr) d / 2 



so 



s(x)\ < E 



7T 



■h+-+3d 



ii! • ••• • id! 



J{n k ,x) 



k=r 



(27T) d / 2 



For brevity denote the outer summand above by hj 1: ___j d (t). Then 



2 \ 2 

dx 



so that 



V,s 2 < 



E 

(ji,- j'd)eJ 



(ji,- ,id)eJ 



£ E 





< 


7PH+-- 




ii!-.. 


• • id! 


vr ii+-- 




ii!-.. 


• • id! 


vr ii+-- 





E i ■ • • • ■ 5 ; 



2 \ 2 

dx 



■id e 
kd 



i{n k ,x) 



(2tt)<V 2 



2 1 



k=r 



ii! • • • • - id! 



7T 



sup 

r<k<s 



E 

(ji,- ,id)eJ 

d / oo /_ sup 



E'° fc 

. fc=r 



9 / sup „ „ \ 2 0i+---+Jd) 

r < k < s 



IK ~ ifelloo) 



Ji! • • • • ■ 3S 



n e 



(^r<T< a K-*fc ll°°) 

J>! 



(Em 2 ) ! 
Em 2 ) 



k=r 



= ( e * rf (r^.ll»*-**H-)_i)(^| ajfc |2)'. 



□ 



Corollary 5.4. Let (nfc)fc g N 6e an enumeration ofL d , and let (tfc)fceN C M d such that 

sup ||n fc - tfclloo = L < oo. 
fceN 

Define ek,fk : M. d ^ C by ek(x) = 7^572 e^ TOfc '^ and 7^x572 e l ^ tk ' x h Then the map T 
L2[— 7r, vr] d — >■ ^2[— vr,7r] d ; defined by Te n = e n — f n , satisfies the following estimate: 
(20) ||T|| < e"" - 1. 
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Proof. Lemma ()5.3p shows that T is uniformly continuous on a dense subset of the ball in 
L 2 (E), so T is bounded on L 2 [— vr, n) . The inequality ([20]) follows immediately. □ 

Corollary 5.5. Let (rik)keN, (^fc)fceN C M. d , and let e^, fk and T be defined as in Corollary 
\5.4\ For each I £ N ; define T/ by Tie/- = e^ — fk for 1 < k < I, and Tie^ = for I < k. If 
limfc-^oo ll^fc — tfc ||oo = 0, then lim^oo T\ = T in the operator norm. In particular, T is a 
compact operator. 

Proof. As 

oo oo I 

(T-2])(^Ja fc e fc ) = ^a fc (e fc - / fc ) - ^a fe (e fc - 

k=l k=l k=l 

oo oo 

= ^2 a k( e k - fk) = T( ^ a kdk)i 
k=l+l k=l+l 

the estimate derived in lemma (|5.3p yields 

oo oo oo 

||(r-T0(^o fe e fc )|| 2 = \\T{ a fc e fc )|| 2 <(e 7r ^>" P +1 |l4!l ° c -l)||^a fe e fc || 2 , 

k=l k=l+l k=l 

so || (T — Ti) || — > as / — > oo. As Ti has finite rank, we deduce that T is compact. □ 
We are ready for the proof of Theorem 15.11 



Proof. Step 1: B is a compact perturbation of the identity map, namely 
(21) B = I + lim (-^ + (PtB- 1 ^)- 1 ). 

I— too 

Since (fk)keN is a Riesz basis for L^y— 7r,7r] d , L* = (I — T) is an onto isomorphism where 
^ e fc = e k ~ fk] so B simplifies to (/ — T) _1 (/ — T*) _1 . We examine 

p- 1 = (I -T*)(I -T) = I + (T*T — T — T*) := I + A, 

where A is a compact operator. If an operator A : H — > H is compact then so is A*, hence 
P/AP/ — > A in the operator norm because 

H^APi-AH < ||PjAPj-PjA|| + ||P,A-A|| < ||APj-A|| + ||P,A-A|| 

= |[P ; A* - A*|| + ||P;A - A|| 0. 

We have 

B- 1 = lim (/ + PiAPj) = lim (/ + P^B' 1 - I)Pi) = lim (/ - P t + PiB^Pi). 

I— too I— too I— too 

Now (PiB~ 1 P[) restricted to the first / rows and columns is the Grammian matrix for the set 
(/li ' ' " > fl) which can be shown (in a straightforward manner) to be linearly independent. 
We conclude that PiB~ 1 Pi is invertible as an / x / matrix. By (P^P" 1 ^) -1 we mean the 
inverse as an I x I matrix and zeroes elsewhere. Observing that the ranges of P[B~ l Pi and 
(P[B~ l P/) _1 are in the kernel of I — Pi, and that the range of / — Pi is in the kernels of 
P;P _1 P/ and (P;P _1 P;) _1 , we easily compute 

{I-Pl + (PiB^Pi)- 1 )- 1 = I-Pi + PiB^Pu 
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B 



-i 



lim (I-Pi + (PiB-'Pi) 

l— >oo 



1 D.\-l\-l 



lim (I-P1 + (PiB^Pi)- 1 ) : = lim B l = I + lim (-Pi + (P^- 1 P^ 1 ) . 



I— >oo 



I— >oo 



I— >oo 



Step 2: We verifiy equation (fT6|) and its convergence properties. Recalling equation (fTTj) . 
we have 



/(*) 



implying 



00 00 
i ~Pl + (PiB-'P^fs/^git -{) = y d E^ ~ B i)fs/xhg(t - -f) 



k=l 



k=l 



fit) 



k=l 
00 



A'' 



fc=l 



Therefore, 



" aJz EK^^^s/aM • -y) 



fc=i 



A 7 II2 



/c=i 



fc=Z+l 



-A7T,A7r] G 



1 

A 5 



^ 1 (»)(-)(E[( b - b ')/5/a]*^ > + e /(y)^ 

fc=l fc=2+l 



-A7T,AT] d 



after taking the inverse Fourier transform. Now 

/ 



/(•) 



1 



Y.KPiB-'Pir'fs/xUi 



k=l 



^ 00 

^ xuHE^-^a]^''" 



A d 
M 

£ 1? 



fc=l 



[-A7r,A7r] d A d 



00 

4 E /(t)^ 



fc=/+l 



[-A7r,A7r] d 



(B - Ws/A 



€a(N) A d 



71 JT 00 j 1 



fc=;+i 



since (/fc(^)) fe is a Riesz basis for Li\— Xtt, \7r] d . Since i?; — > -B as Z — > 00 and (/(x))jk e 
^2(N), the last two terms in the inequality above tend to zero, which proves the required 
result. 
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Finally, to compute (PiB~ l Pi) nm , recall that B^ 1 = (I — T*)(I — T). Proceeding in a 
manner similar to the proof of equation fjlQf) , we obtain 

^mn (LL e n ,e m ) (L e n , L e m ) {fni fm) 

= sinc7r(t nj i - t m> i) ■ ... ■ sinc7r(i nid - t m ^ d ). 
The entries of PiB~ 1 Pi agree with those of B~ l when 1 < n, m < I. □ 

One generalization of Kadec's 1/4 theorem given by Pak and Shin in |12| (which is actually 
a special case of Avdonin's theorem) is: 

Theorem 5.6. Let (i n ) n gz C M. be a sequence of distinct points such that 

limsup |n — t n \ = L < —. 

Then the sequence of functions (fk)k&, defined by fk{x) = ^j^ e%tkX ' > *s a Riesz basis for 
L 2 [-TT,ir]. 

Theorem 15.61 shows that in the univariate case of Theorem 15.11 the restriction that (ifc)fceN 
is a Riesz basis for L2\— vr,7r] can be dropped. The following example shows that the mul- 
tivariate case is very different 

Let ( 

e n)n be an orthonormal basis for a Hilbert space H. Let f\ £ H with ||/i|| — 1, then 
(/i> 62, e3, • • • ) is a Riesz basis for H iff (fi, e\) ^ 0. Verifying that the map T, given by 
efc >■ ct for k > 1 and e\ i— > /i, is a continuous bijection is routine, so T is an isomorphism 
via the Open Mapping Theorem. In the language of Theorem I5.R (/i, e2, e%, • • • ) is a Riesz 
basis for L,2[— n, n] iff 

^ sinc(7rti 5 i) • . . . • sinc(7rt 1 
that is, iff ti G (R \ {±1, ±2, • • • 



6. Generalizations of Kadec's 1/4 Theorem 
Corollary 15.41 yields the following generalization of Kadec's Theorem in d dimensions. 

Corollary 6.1. Let (n^)^^ be an enumeration ofL d , and let {tk)k&i C M d smc/i £/ia£ 

, , „ . ln(2) 

(22) sup ||n fc - tfclloo = £ < — v-> 

fceN 

T/ien i/ie sequence (fk)keN defined by fk{ x ) = ^ 2 ^ d y 2 e^ :r '* fc ^ is a Riesz basis for L2[— it, Tr] d . 

The proof is immediate. Note that equation (|20p implies that the map T given in Corol- 
lary [53] has norm less than 1. We conclude that the map (I — T)e^ = fk is invertible by 
considering its Neumann series. 

The proof of Corollary (|5.4p and Corollary (|6,ip are straightforward generalizations of the 
univariate result proved by Duffin and Eachus [8j. Kadec improved the value of the constant 
in the inequality (|22p (for d = 1) from ^jp- to the optimal value of 1/4; this is his celebrated 
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"1/4 theorem" [10]. 

Kadec's method of proof is to expand e tSx with respect to the orthogonal basis 

{1, cos (nx), sin (n - ^)x}„ eN 

for L^\— 7T, tt], and use this expansion to estimate the norm of T. In the proof of Corol- 
lary (|5.4p and Corollary (|6.ip we simply used a Taylor series. Unlike the estimates in 
Kadec's Theorem, the estimate in equation ([20]) can be used for any sequence (tk)ken C R rf 
such that sup fcgN \\nk — tfclloo = L < oo, not only those for which the exponentials (e ltnX ) n 
form a Riesz basis. An impressive generalization of Kadec's 1/4 theorem when d = 1 is 
Avdonin's "1/4 in the mean" theorem, [I]. 

Sun and Zhou (see [13] second half of Theorem 1.3) refined Kadec's argument to obtain a 
partial generalization of his result in higher dimensions: 



Theorem 6.2. Let (a n ) ngZ d C R such that 

0<L<\, 

/ . sin7rL\ rf /'simrL\ d 

Ud{L) := I 1 — cos ttL + sm7rL H — I — — , and 

V ttL J V ttL J 

||o n — n||oo < L, n G Z rf . 

If Dd{L) < 1, f/ien ( e ^° n ) * s a ^ esz basis for L2\— 7r, 7r] d ttfit/t frame bounds (1 — 

Dd{L)) 2 and {l + D d {L))\ 

In the one-dimensional case, Kadec's theorem is recovered exactly from Theorem 16.21 When 
d > 1, the value satisfying < < 1/4 and Dd(xd) = 1 is an upper bound for any 
value of L satisfying < L < 1/4 and Dd(L) < 1. The value of x,i is not readily apparent, 
whereas the constant in Corollary 16. II is ^j. A relationship between this number and Xd is 
given in the following theorem (whose proof is omitted). 

Theorem 6.3. Let x^, be the unique number satisfying < x^ < 1/4 and D^(x^) = 1. Then 

x,- x ^ 
lim — ?£. = 1 

i^o (In 2)2 
12-kcC 2 

Thus, for sufficiently large d, Theorem 16.21 and Corollary 16.11 are essentially the same. 



7. A METHOD OF APPROXIMATION OF BIORTHOGONAL FUNCTIONS AND A RECOVERY OF 

A THEOREM OF LEVINSON 

In this section we apply the techniques developed in the previously to approximate the 
biorthogonal functions to Riesz bases (—k=e ltn ^) for which the preframe operator is small 

^ V Z7T ' 

perturbation of the identity. This is the content of Theorem 7.1. A well known theorem of 
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Levinson (see pages 47-67]), follows as a corollary to Theorem 7.1. 



Definition 7.1. A Kadec sequence is a sequence (t n )nez of real numbers satisfying 

sup \t n — n\ = D < 1/4. 

Theorem 7.2. Lei (£ n )nez C R fre ci sequence (with t n ^ /or n ^ 0) such that (/n)n — 
(-y==e jtn (''') n is a Riesz basis for L2[—7t,tt], and let (e n ) n be the standard exponential or- 

thonormal basis for L2[— 7r, 7r]. // £/ie map L given by Le n = f n satisfies the estimate 
\\I — L\\ < 1, then the biorthogonals G n of -i=J 7 (/ n )(-) = sinc(-7r(- — t n )) in PW[_ n ^ are 



(23) ^ i-w ' " £Z - 

w/iere 

oo 

(24) Hw-d-^nCi -j-)(i -^). 

n=l 

Definition 7.3. Let (i n )nez Clbea sequence such that (f n )n = (^|^ ei *™^ ) n ^ s a Ri esz 
basis for Li\— tt, 7r]. If / > 0, the /-truncated sequence (ii, n )nez is defined by fy n = i„, if 
1 77, | < / and ti n = n otherwise. Define // n = ^=e* ii ' ,I ^'- ) for n G Z, I > 0. 



/2tt 

Let Pi : Ij2[—tt, it] — > L2[—tt, ir] be the orthogonal projection onto span{e_;, . . . , e{\. 



Proposition 7.4. Let (t n ) n & CM be a sequence such that (f n ) n (defined above) is a Riesz 
basis for L2[— tt, tt]. If (e n ) n is the standard exponential orthonormal basis for L2[— vr, ir] and 
the map L (defined above) satisfies the estimate \\I — L\\ = 5 < 1, then the following are true: 

1) For I > , the sequence (fi, n )n i> s a Riesz basis for L2[—tt,it]. 

2) For I > 0, the map Li defined by Lie n = fi jTl satisfies \\Lf 1 \\ < j^r. 

Proof If (cn) n € £ 2 (Z), then 

(I - Li){y y c n e^) = ^c n (e„ - L ; e„) = ^(e„ - /„) = (/ - L)P;(^c n e n ), 

n n \n\<l n 

so that 

(25) (/ - L() = (I - L)P h 

From this, ||J — Li\\ < 5, which implies 1) and 2). □ 

Define the biorthogonal functions of (fi <n )n to De (/* n ,)n- Passing to the Fourier transform, 
we have -^T(f l)n )(t) = sinc(7r(i and Gi >n (t) := -^=J r (/* n )(t). Define the biorthog- 
onal functions of (f n )n similarly. 



so that for all m, / z * m = (Xf) 1 e m . Similarly, = (L*) 1 e m . We have 
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Lemma 7.5. If (t n ) n C M satisfies the hypotheses of proposition \7J\ then 

lim Gi n = G n 

in PW { _ n ^]- 
Proof. Note that 

<5nm = \fl,n>fl, m ) = \^l e ni fl, m ) = ( e n, L*l fl iT n) 

i,m = m 

Now equation (|2"5j) implies L — Li = (L — /)(/ — Pj), so that 

/i% - & = (^TV - Pi)(L* ~ I){L*Y l e m . 
Applying proposition 17.41 yields 

11//% " fraW < J^gW ~ Pl)(L* ~ I)(L*r l e m l 

which for fixed m goes to as I — > oo. We conclude linn^oo /* = which, upon passing 
to the Fourier transform, yields lim^oo G/, m = G m . □ 

Proof of Theorem 17.21 

We see that 5 nm = (Gi >m , Si >n ), where Si >n (t) = sinc(7r(t — t n )) when \n\ < I and Si >n (i) = 
sinc(7r(t — n)) when \m\ > I. Without loss of generality, let \m\ < I. Equation ([1]) implies 
that G[ m (k) = when \k\ > I. By the WKS theorem we have 

k=l k=l 

Gi, m (t) = Gi, m (k) S mc(7r(t - k)) = ( £ ^^)sinc(^) 

k=-l k=-l 

= — ; smc(7rt), 

nU(* -*)(-*-*) 

where W[ is a polynomial of degree at most 21. Noting that 

00 f 2 1 1 f 2 

sinc(vrt) = J] (1 " tj) and l[(k - t)(-k - t) = (I!) 2 ]J (l - ^) , 
k=l ' k=l k=l 

we have 

r m (-i)W*) n h *\ 

V '' k=l+l 

Again by equation ([1]), 5 nm = Gi^ m (t n ) when \n\ < I so that 

(— 1V °° 
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This determines the zeroes of W[ . We deduce that 

,„ m _ ciUt=i(t-t k )(t-t- k ) 
^ > ~ TZT+ 

for some constant q. Absorbing constants, we have 
where 

k=l k k l+l 

Now = Hi(t m ), so G|, m (t) = Q g ' (t ilg (tm) - Taking limits, q = (H; y (tm) - This yields 

Gl, m (t) - 



Define 

^(t> = - i0 ) n (! - f ) c 1 - A)- 

k=l k ~ k 

Basic complex analysis shows that H is entire, and Hi H and Hi — > H' uniformly on 
compact subsets of C. Furthermore, H'(tk) 7^ for all k, since each t k is a zero of H of 
multiplicity one. Together we have 

By the foregoing lemma, G^ m — > G m . Observing that convergence in PW t j implies 
pointwise convergence yields the desired result. 

Levinson proved a version of Theorem l7.2l in the case where {t n ) n& z is a Kadec sequence. His 
original proof is found in [Til pages 47-67]). We recall that if (f n )n is a Riesz basis arising 
from a Kadec sequence, then the preframe operator L satisfies ||J — L\\ < 1. Levinson's 
theorem is then recovered from Theorem 17.21 
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